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Satellite formation missions require the precise determination of both the position and 
attitude of multiple vehicles to achieve the desired objectives. In order to support the 
mission requirements for these applications, it is necessary to develop techniques for rep- 
resenting and controlling the attitude of formations of vehicles. A generalized method for 
representing the attitude of a formation of vehicles has been developed. The representation 
may be applied to both absolute and relative formation attitude control problems. The 
technique is able to accommodate formations of arbitrarily large number of vehicles. 

To demonstrate the formation attitude problem, the method is applied to the attitude 
determination of a simple leader-follower along-track orbit formation. A multiplicative 
extended Kalman filter is employed to estimate vehicle attitude. In a simulation study 
using GPS receivers as the attitude sensors, the relative attitude between vehicles in the 
formation is determined 3 times more accurately than the absolute attitude. 

I. Introduction 

T here are many satellite formation flying missions in which the coordinated control of the formation’s 
attitude is proposed . * 1 For example, it is necessary to reposition and reorient the satellites in a variable 
baseline interferometer application to image different targets. The size and cost of the individual satellites 
in the formation are commonly constrained to be as small as possible for the given mission requirements. 

For Earth orbiting formations, it is efficient to use Global Positioning System (GPS) receivers as spacecraft 
sensors. The receivers can provide autonomous on-orbit position determination and time synchronization 
between satellites. Current GPS receiver technology conserves power and is well-suited to miniaturization. 
Dynamic filters can be used to provide estimates even when measurements are temporarily unavailable. In 
this manner, the useful range of spaceborne GPS receivers can be extended to altitudes above the GPS 
constellation, increasing the diversity of applications that can benefit from GPS devices. 

It is also possible for GPS receivers to determine a vehicle’s attitude. There are now several different 
methods by which the attitude information may be obtained. The most common method is a local interfer- 
ometric one using multiple antennas on a single vehicle . 2 This technique requires the resolution of a carrier 
phase cycle ambiguity by external or other means . 3,4 Other methods employ antenna gain pattern mapping 
with signal to noise ratio measurements . 5 It is possible to use this approach with a single GPS antenna and 
another sensor, such as a magnetometer, to determine three-axis attitude . 6 

Since a GPS sensor will almost certainly be present on an Earth orbiting satellite for position and time 
determination, it would be advantageous to employ that device for attitude determination if it can reliably 
meet the mission attitude requirements. These requirements could be given in either absolute or relative 
attitude determination accuracy, depending on the mission. 

It is well known that the relative position solution between two GPS receivers located in near proximity 
can be obtained more accurately than the absolute solution . 7,8 The greater accuracy is possible because 
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the most significant error sources are correlated between the different GPS measurements. When the range 
measurements are differenced, a more accurate relative solution is directly obtained. 

For vehicle formations using GPS receivers, it might be expected that attitude measurement error can- 
cellation would also occur between similarly designed vehicles operating in proximity with each other. If the 
measurement errors are highly correlated, relative attitude determination may be performed more accurately 
than absolute attitude determination. 

Survey of Prior Work and Outline of Topics 

Although formation orbit determination and control has been extensively studied, the formation attitude 
determination problem is somewhat less established. Scharf et al have compiled a survey of formation flying 
guidance and control applications with more than 100 references. 9 The papers that address formation atti- 
tude tend to focus primarily on control methods. 10-14 This paper approaches the problem from a different 
perspective: the formation attitude control is assumed and the sensor model is developed. First, a generalized 
framework is presented in which the formation’s dynamics may be expressed. Then a quaternion attitude 
estimation algorithm is developed using GPS carrier phase measurements. The generalized framework and 
estimation algorithm are demonstrated using a Low Earth orbit two vehicle leader-follower formation ex- 
ample. The effect of having GPS correlated measurement errors is evident in the resulting improvement in 
relative attitude accuracy that is achieved. 


II. Reference Frame Definitions for Vehicle Formations 


To describe the dynamic state of a formation of vehicles, some conventions must be established in the 
definition and notation of frames of reference. The approach used in this paper is similar to one that was 
previously introduced by Xing. 15 The main feature of this approach is that it is generalized and scales easily 


with the number of vehicles in the formation. 

The fundamental reference frame shown in fig- 
ure 1 is an inertial reference specified as Fj. The for- 
mation target reference frame is given as Ft, which 
may be non-inertial. This frame defines a local ref- 
erence for the entire formation. It does not have to 
be occupied by an actual satellite. Fn is the target 
reference frame for the Pth satellite in the formation. 
This is considered to be the satellite’s desired posi- 
tion and orientation relative to the formation target 
reference frame. Finally, the Pth satellite’s actual 
position and orientation is given by the Fbi body 
reference frame, which is attached to the Pth vehi- 
cle’s center of mass. In a regulator problem, the goal 
is to align Fsi with Fn within some acceptable tol- 
erance. A general vector v may be coordinated with 
respect to the Pth local body frame, Fbi- 

In the notation of figure 1 ( the translational ve- 
locity of the Pth body frame relative to the Pth tar- 
get frame is simply given by 



Tl d 

dt 


( A Pj 


(1) 


where A p. is the vector from the Pth target frame to the Pth body frame. The velocity of the body 
frame to the formation target frame can be expressed by adding translational and rotational terms 


T _d 

dt 


(BbJ 




( 2 ) 


In this case, uj TTi is the angular velocity of the Pth target frame (Fn) relative to the formation target 
frame (Ft). The body frame’s (Fsi s) velocity can be similarly expressed relative to the inertial frame (Fj) 
with the addition of two more terms 
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( 3 ) 


r _d 

dt 


(r B i) 


Ti, 


x £T, + ^( A a)+H - * A £ 


The translational accelerations may be likewise derived. With the translational portion of the problem 
thus stated, it will no longer be considered in this paper, although its presence in the coupled formation 
control problem is assumed. 

The attitude of a reference frame may be related to another by a unit quaternion, which is based upon 
Euler’s theorem 


q = 


q 


esin (0/2) 

. 94 . 


cos (0/2) 


( 4 ) 


where e is the 3x1 axis unit vector and 0 is the angle of rotation that relates the two frames. 
Employing the definition of quaternion multiplication as stated by Shuster, 16 


qA®q B = 

the quaternion conjugate, 
and the vector quaternion, 


(q Ai q R + q Bi q A -g A x g R ) 


B ’ 


q Ai q B 4 -g A -g B 


( 5 ) 


q = 


1 T 


<74 


_ 5 


( 6 ) 
( 7 ) 

vn = qri ® v B i <X> qri (8) 

which takes a vector quaternion expressed in one frame and transforms it to another frame. Eq. (8) is 
equivalent to the well known direction cosine transformation matrix 


the quaternion rotation operation is defined 


Bi , 


1 T 


-Bi 


—Ti 


— fjBi 

— ° Ti V-Bi 


and in fact may be derived from q®\ as follows, 17 

C(q) = ( q\ - M 2 ) • hx 3 - 2 q 4 [gx] + 2 q q T 
where I^ x ^, is the 3x3 identity matrix and [gx] is the cross product matrix 


(9) 

( 10 ) 


0 -g3 q2 

[gx] = q 3 0 — gi (11) 

Q '2 gi 0 

Successive rotations are defined using the quaternion rotation operator listed in eq. (8). For example, 


Bi 

VT = q T ® VBi 


, g ^ = qdp 1 ( 


q?t 1 


' VBi ® qrl ® <7t* 


(12) 

It is evident that the quaternion g|?| is a representation of the attitude of the i’th satellite local body 
reference frame (F B i) relative to the i’th satellite target reference frame (Fn). If the target reference frame 
tracks the desired satellite state, then the generalized vehicle attitude tracking problem can be given as 


„ Bi 


In ’ ^ 


TiBi 


o 

(q^Ti 4 =i 


as t 


(13) 


There is a sign ambiguity on the scalar part of the quaternion in equation (13). An attitude performance 
vector for the i’th satellite may be established as 


Zi(i) = 


1 T 


M - 1 


(14) 
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Formulation of Optimal Attitude Controller 

Using the performance vector defined in eq. (14), a regulator performance index may be specified as (for 
example) 


Mt) = 



(zf QiZi + uj RiUj)dt 


(15) 


which incorporates u, ( as the actuator control effort, and Qi and Ri are gain weighting matrices. z t may 
also include translational states if a coupled translational and rotational controller is required. Although not 
shown, Ji may include operational constraints using Lagrange multipliers and the calculus of variations. 
The locally optimal controller C,; at the i’th satellite may be specified as 


C i : min Ji(t) (16) 

t — too 

Over the entire formation, the globally optimal controller C is given by 

C : min Ji{t) (17) 

i 

This completes the statement of the optimal control problem using the general framework presented. 
From this point forward, the paper focuses on the attitude determination problem, and the existence of a 
sufficiently optimal controller is presumed. 


III. Attitude Estimation 

The method of attitude representation chosen for this algorithm is the quaternion, as given in eq. (4). 
The quaternion is generally preferred for spacecraft applications because of its lack of singularities and its 
computational efficiency. However, there are only three degrees of freedom in the rotation group and so 
precautions must be taken to prevent the covariance matrix from becoming singular over time. 

Following the method employed by Markley, 17 the true attitude is represented as a small error deviation 
(< 5q) from a reference attitude quaternion ( q re f ), as expressed by the quaternion product 

q(t) = 6q{a{t)) 0 q re f(t) (18) 

Note that the deviation vector a(t) has only three degrees of freedom and the components of a are called 
Modified Rodrigues Parameters 18 


a = 4e tan (0/4) 


(19) 


The factor of 4 causes the magnitude of |a| « 0 for small rotations. The second order approximation to 
the error quaternion is given by 17 


Sq(a) 


a/2 

1 — | o | 2 /8 


( 20 ) 


This method, known as the multiplicative extended Kalman filter (MEKF), avoids numerical stability 
problems that can occur from the normalization of the quaternion q. 


Time Update 

Attitude state propagation is performed using the kinematic quaternion relation 


1 

Qref — ^ 


—ref 

0 


0 qref 


( 21 ) 


If the vehicle’s attitude motion approximates that of a rigid body in an inertial reference frame, Euler’s 
equations may be used to propagate u_ re j 


ikref ~ J 1 \(*J—ref) X ^-ref + M_D + Me] 


(22) 
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where it is assumed that the body frame (Fg,) defines the principal axes of the vehicle. J is the 3 .t 3 
diagonal matrix of principal moments of inertia, and M D and M_ c are disturbance and control moments, 
respectively. 

Covariance propagation is performed according to the differential equation 


P = FP + PF t + GQG t 

The linear approximation to the time derivative of a is 17 

a = /(a, t) ~ —u re f x a + 77(f) 

and 77 (f) is a zero-mean white noise process. 

The matrices listed in eq. (23) are then given as 

O p 

F W =da = - [—re/ x] 

df 

G(t) = — = I 3x 3 

and 


E \jl(t)rj r ( t)\ = Q(t)6(t - t) 


(23) 

(24) 

(25) 

(26) 

(27) 


Measurement Update 

Up to this point the filter design has been independent of the sensing technology that is needed to determine 
the attitude of the vehicle. If the vehicle is in Earth orbit, it is possible to sense its attitude with a GPS 
receiver. 

The method of attitude determination employed 
in this study is GPS interferometry using differen- 
tial carrier phase measurements. It is assumed that 
the cycle ambiguities have been resolved, therefore 
the carrier phase measurements may be converted 
directly to range A r measurements between anten- 
nas. The geometry of this measurement is shown 
in figure 2. The carrier phase measurement A <j>ij is 
made between two antennas with baseline vector b-j 
tracking a GPS satellite with line of sight vector l t 


A(f>ij — A i'ij — kij + Pj + Vij — Zj ■ bj — kij + Pj + Vij 

. (28) 

kij is the known single difference cycle ambiguity, 
Pj is a bias term, and i'ij is measurement noise. The 
bias term is removed by differencing measurements 
along the same baseline j between GPS satellites i and k 

U Ao'),j/,: A pij A (fikj — bj ' (b f p, ) ( k-jj Zi/.y j -(- (i'ij 7'7/;:y ) (20 j 

Since the integers k are known, they may be incorporated into the carrier phase measurements and the 
measurement equation may be rewritten as 

V = bj- (h - Ik) + (vij ~ Vkj) = h + u (30) 

In order to evaluate the dot product in eq. (30), the vectors must be coordinated in the same reference 
frame. The line of sight vectors (/) are known in the external reference frame, but the antenna baseline 
vectors (b) are given in the body frame, so they must be transformed to evaluate the dot product 



Figure 2. Carrier Phase Interferometry. 
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h=[C(q)-b j ] T ■ ({,-y (31) 

The quaternion q is now substituted using equations (18), (10), and (20) to give the linear approximation 

C(q) = C(5q <g> q re f) = C(6q)C(q ref ) « (I 3x3 - [ax]) C(q re f) (32) 

h may now be written in terms of the attitude deviation vector, a 

Ha) ~ [(-^ 3 x 3 - [«x]) C{q ref ) ■ b j \ T - (Z 4 - l_ k ) (33) 

which may be rearranged as 

h(a) Kh+ {li - l_ k ) T ■ ( C{q ref ) ■ bj) x a (34) 

where h is the expected measurement at the attitude given by q re f • The measurement sensitivity matrix 
is 

n T = Qa = - l*) x ] • ( C ^ef) ■ b_ 0 ) (35) 

The Kalman Gain matrix is computed as 

K = P{-)H T [HP(-)H T + R]- 1 (36) 

where R is the covariance of the zero mean measurement white noise, E \v(t)v T {t)\ = R(t)S(t — r). The 
error state update is then completed 


a = K(y - h) 

(37) 

Qref(+) = Sq(a) <g> q re f{—) 
The error covariance is updated according to 

(38) 

P(+) = (/ 3 * 3 - KH) ■ P(-) 

(39) 


This completes the derivation of the vehicle attitude estimator. The filter will be adjusted for relative 
attitude observations in a later section. 

IV. Relative Attitude Simulation 

For the purposes of demonstrating the relative attitude concept using GPS receivers, a 6 N degree of 
freedom simulation was developed, where N is the number of vehicles in the formation. In the example used 
in this paper, N = 2, however the results can be generalized to N > 2. 

The translational dynamics are specified as a simple leader-follower arrangement with 2 satellites in the 
same orbit whose elements are specified by table 1 The anomaly of the follower satellite was adjusted to 
lag the leader satellite by a constant distance, R = 100 km. Since the focus of this study is on the attitude 
estimator performance, no perturbations are considered in this basic simulation. 

Table 1. Orbit elements for leader satellite 
at t = 0.* 

Element Value Element Value 

h 800.0 km Q 0 

i 89.9 deg u 0 

e 0 r'o 0 

* Follower satellite lags leader by 100 km. 
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The formation target reference frame (Ft) is assigned to the leader satellite with the 
z-axis pointed along the zenith vector, the j- axis is rotated 90 degrees in the direction 
of motion, and ixj = k. The frame rotates with the leader satellite’s center of mass, 
and therefore defines a Local Vertical Local Horizontal (LVLH) reference for this circular 
orbit. The follower satellite is assigned a similarly defined LVLH target reference frame 
at its location (i 7 ^)- The attitude of each vehicle’s body axes (Fbi and Fbi) are given 
relative to their respective target reference frame (Ft and Fr-i)- Since the satellites are 
not in the same position, Ft ^ Ft 2 - 

The rotational dynamics are specified as follows. A 50 kg, axisymmetric satellite is 
proposed as shown in figure 3 The satellite is a twin-boom configuration with each boom 
extending 3.05 m (10 ft) in length. The main body of the satellite is a 40 kg homogeneous 
cylinder with radius 0.23 m and height 0.38 m. The booms are modelled as lumped masses 
of 2.5 kg at the center of each span. Each boom tip is also modelled as a homogeneous 
cylinder, with mass 2.5 kg, radius 0.115 m, and height 0.19 m. The entire structure is 
assumed to be rigid and symmetric. The principal axes of the structure define the local 
body frame (F B ). 

Using this simple mass model, the inertia properties are calculated as J x = 1.09 kg m 2 , and J y = J z 
kg m 2 . The non-dimensional inertia ratios are 


25 kg 


Figure 3. 
Gradient 
lite. 


25 kg 

Gravity 

Satel- 


= 71.4 


kpitch — {Jy Jx) / Jz — 0.985 


(40) 


kroll = (Jz ~ Jx) / Jy = 0.985 


(41) 


v yaw 


— (Jz — Jy) / Jx — 0 


(42) 


According to linearized gravity gradient dynamic analysis, the vehicle attitude motion is represented with 
respect to the LVLH frame according to the following differential equations 19 


Yaw: 

Roll: 


1 

u 

^ yaw 

l 


^ roll L 


Ipx - nipy 
i’y + nipx 


nipy + n 2 ip x 


nip x - An 2 ip v 


0 + (M cx + Mdx) 

= 0 + (M cy + Mdy) 


(43) 

(44) 


Pitch : 


k 


pitch - 




[3n 2 ip z 


0 + (M cz + Mdz) 


(45) 


n is the angular rate, which is 0.595 rev/hr for this case. Control and disturbance moments could 
be applied to the right hand side of eqs. 43- 45 as shown, but will be assumed 0 for this analysis. The 
homogenous solution to the linearized gravity gradient equations of motion for small deflections is a pair of 
oscillators. For this configuration, the natural frequencies are 


u plt ch = n(3k pltch ) 0 - 5 = (0.595)(3 • 0.985) 0 ' 5 = 1.02 rev/hr (46) 


u> ry = n( 1 + 3 kroll ) 0 - 5 = (0.595)(1 + 3(0.985))°' 5 = 1.18 rev/hr (47) 

The unforced response of the leader satellite is simulated with an initial condition offset of 30 degrees 
in pitch and roll relative to the target reference frame (q^ 1 )- This response is shown in the left hand side 
of figure 4 on the next page. The follower satellite undergoes a very similar motion which is offset in phase 
angle by the difference in anomaly. The true relative attitude between the two vehicles is given by (see eq. 12 
on page 3) 


9bi = 9m 


' 9 ? : 2 = 1 


®qf 2 


(48) 


The relative attitude motion between the leader and the follower is shown as a 3 — 2 — 1 Euler angle 
sequence in the right hand of figure 4 on the next page. The Euler angles are extracted from the quaternion 
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40 


Leader Satellite, Gravity Gradient Attitude Simulation 


Leader-Follower Relative Attitude, Expressed in F B1 Frame 




(a) Leader Simulated Attitude. (b) Leader-Follower Relative Attitude. 


Figure 4. Simulated leader absolute attitude (q^ 1 ) and leader- follower relative attitude 


q B \. Since the satellites are only 100 km apart, the relative attitude motion between the leader and follower 
satellites is small. The pitch bias is due to the difference in local vertical direction between the two satellites. 

It is possible that the initial conditions of the gravity gradient motion between the satellites could be 
offset by different amplitudes and phase angles. In this case, there would almost certainly be larger relative 
motion between the two vehicles. However, due to the fact that the satellites have the same mass properties 
and orbits, it is expected that over time they will settle into gravity gradient motion of similar amplitude 
that is phase shifted by the value of the anomaly. This case is the one that is simulated. 


V. Relative Attitude Results Using GPS Receivers 


To simulate the performance of a GPS receiver in this application, the attitude estimator in section 
III processed simulated GPS measurements for each satellite undergoing the dynamics in section IV. The 
estimator is given simulated carrier phase measurements for all GPS signals that are within a 60 degree cone 
angle for each satellite’s upward-facing axis of symmetry (the i Bl and i B2 directions). The carrier phase 
measurements are made from GPS antennas that lay along the vertices of a 1 m cube that is aligned with 
the satellite’s body axes (3 total baselines per satellite). All antenna boresight directions are pointed along 
the vehicle axis of symmetry. Signal blockage due to structure is not considered, but carrier phase multipath 
is simulated in the form of correlated measurement noise. 

To compute the relative attitude, an estimator is run independently for each satellite, producing the 
absolute attitude estimates, qf 1 and qf 2 . The relative attitude estimate is then computed (compare to 
eq. 48 on the preceding page) 


9b\ = 9m 


' qf 2 = qf 1 ® qf 2 


(49) 


The resultant attitude error is then computed as follows. For the absolute attitude error of the leader 
satellite, the leader absolute error quaternion is given by 


— r , 1 Z B1 — Z B1 

qabserr — 9b 1® 9l — 9 1 ® 9l 

The relative attitude error of the follower satellite to the leader satellite is likewise specified 


(50) 


Qrelerr = qf\ ® 9m = 9 b1 ® 9b1 ( 51 ) 

In each case, the error is reported as a 3 — 2 — 1 Euler angle sequence, where the error Euler angles are 
extracted from the error quaternions q a bserr and q re ierr- Performance statistics are computed using the error 
Euler angles. 
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(a) Leader Attitude Error. 


(b) Leader-Follower Relative Attitude Error. 


Figure 5. Noise free absolute attitude error (q a b serr ) and leader-follower relative attitude error ( f/ ./ ..... ) . 


Noise Free Simulation 

Before introducing the effects of measurement noise, it is instructive to observe the estimator’s performance 
in the noise free case. This may be considered as a best case bound on the performance for this algorithm. 

The result, shown in figure 5. demonstrates the main thesis of this paper. The three-axis absolute attitude 
performance of 0.181 degrees root mean squared (RMS) is typical for this sensor under these dynamic 
conditions. There is a slight bias of absolute attitude in pitch. The estimator is limited by geometrical 
considerations associated with the number and direction of available measurements at each measurement 
epoch. Also, although the measurements are perfect, the estimator is not designed for this case. The absolute 
attitude performance could be improved by “tuning” the estimator for this special case, but this is not done 
since the actual application includes measurement noise. 

The relative attitude accuracy is strikingly better at 0.006 degrees RMS. Although the absolute attitude 
estimators for both satellites have errors associated with dynamics and measurements, these effects largely 
cancel out and leave a much more accurate relative state. The jumps that occur in figure 5o are associated 
with rising or setting GPS satellites. The jumps could be removed by adding logic to handle these cases. 

Uncorrelated Measurement Noise Simulation 

Every sensor is subject to measurement noise and GPS receivers are no exception. In the case of differential 
carrier phase measurements, the predominant source of noise is signal multipath due to reflections from 
nearby surfaces. Although deterministic in nature, multipath is difficult to model due to the high number of 
possible reflections. The stochastic signature of this noise is that of a time correlated random variable with 
a characteristic standard deviation (cr„) and time constant (t). 20 

In satellite applications, multipath can be a severe source of measurement noise. There are usually a 
large number of highly reflective surfaces in the vicinity of the GPS antenna. Multipath can be minimized by 
locating the antennas carefully to reduce reflections and by utilizing hardware practices such as the addition 
of antenna choke rings. Based on the results of previous orbit experiments, a standard deviation a n = 5 cm 
is considered to be typical for a Low Earth orbit satellite attitude application. 21 

In the worst case scenario of totally uncorrelated measurement noise (correlation coefficient p = 0), there 
is no benefit to differencing attitude solutions between satellites. In fact, relative attitude solution error is 
increased due to the independence of the noise sources. This case is shown in figure 6 on the next page. The 
introduction of measurement noise has increased the absolute attitude error to 0.295 degrees RMS, but the 
most significant development is that the relative attitude error has increased to 0.327 degrees RMS. 
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Leader Absolute Attitude Error, 5 cm uncorrelated measurement noise Leader-Follower Relative Attitude Error, 5 cm uncorrelated measurement noise 




Time (hours) Time (hours) 

(a) Leader Attitude Error. (b) Leader-Follower Relative Attitude Error. 

Figure 6. Uncorrelated noise absolute attitude error ( q a bserr ) and leader-follower relative attitude error ( q r elerr )■ 


Correlated Measurement Noise Simulation 

It is hypothesized that for similar satellite designs that are operating in similar orbits and attitudes, error 
sources that are geometrically dependent such as multipath will be correlated to some extent. In order to 
investigate the effect of correlated measurement noise on the relative attitude solution accuracy, the carrier 
phase measurements are simulated with varying levels of correlated noise between the leader and follower 
satellites. As the correlated noise is added, the total noise is adjusted so that a n = 5 cm. Since the total 
noise level is the same in a probabilistic sense, it is possible to see the benefit of attitude differencing when 
correlated noise sources are present. 

Figure 7 on the following page shows that this is the case. This plot shows the attitude performance 
when the multipath noise is highly correlated (p = 0.9) between satellites. The absolute attitude accuracy 
is approximately unchanged at 0.279 degrees RMS, but the relative attitude accuracy is improved at 0.084 
degrees RMS. 

These results are collected in table 2. An intermediate case for p = 0.5 is also listed. As expected, the 
relative attitude performance improves with the amount of multipath correlation. 

VI. Conclusion 

This paper examined how GPS receivers may be employed as attitude sensors in Earth orbiting satellite 
formations. GPS receivers enjoy advantages as standalone navigation and attitude sensors when the size and 
cost of the vehicles are tightly constrained. It may be possible to meet the mission requirements in some 
cases with a minimal hardware complement of a single GPS receiver and possibly a simple additional sensor 
such as a magnetometer or sun sensor. 

The mission presented here is hypothetical, but it is considered to be representative of several Earth 


Table 2. Performance Summary for Attitude Simulation. 


Case 

Comment 

On 

(cm) 

P 

®abserr 

(deg RMS) 

Qrelerr 

(deg RMS) 

1 

Noise Free 

0 

- 

0.181 

0.006 

2 

Uncorrelated 

5 

0 

0.295 

0.327 

3 

Somewhat Correlated 

5 

0.5 

0.291 

0.210 

4 

Highly Correlated 

5 

0.9 

0.279 

0.084 
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Leader Absolute Attitude Error, 5 cm correlated measurement noise, p=0.9 


Leader-Follower Relative Attitude Error, 5 cm correlated measurement noise, p=0.9 



Time (hours) 



(a) Leader Attitude Error. 


(b) Leader-Follower Relative Attitude Error. 


Figure 7. Correlated noise (p = 0.9) absolute attitude error ( q a bserr ) and leader- follower relative attitude error 

( Qrel err') • 


orbiting formation applications. The design of a GPS attitude estimator is codified using a multiplicative 
extended Kalman filter (MEKF) method. The typical performance of such a system is demonstrated in 
the presence of measurement noise. When the measurement noise is correlated among the satellites in the 
formation (correlation coefficient p > 0.5), improvements in the relative attitude accuracy can be obtained. 
For the case considered in this study, highly correlated GPS measurements (p = 0.9) produced more than a 
threefold improvement in relative attitude accuracy. 

Although GPS measurements will never be fully correlated, it may be possible to create a high level of 
correlation by using the same satellite design and controlling the vehicle attitude so that the geometry of 
the multipath reflections is similar for each satellite. Attitude control is not addressed in this paper, but a 
general framework for evaluating the optimality of a formation attitude control system is given. 

This work represents an early step in the design of formation attitude estimation systems. In this 
study, the GPS measurements were processed separately for each satellite and the attitude solutions were 
differenced. Although good results were obtained, additional work could improve the performance of the 
relative estimation algorithm by centrally processing all of the available measurements. This would allow 
the estimator to take maximum advantage of all the correlations in the measurements. It would also be 
valuable to demonstrate the performance of this system using a real-time hardware-in-the-loop test facility. 
Experience has shown that hardware tests can expose implementation issues that are not apparent in a 
simulation study. These additional investigations are left as future work. 
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